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Abstract 

We compute the 2-point (equal-time) electronic Green function 
in a Tomonaga-Luttinger system with long range electron-electron 
interactions. We obtain an analytical expression for a "super long- 
range" potential of the form V(x) = , As a consistency check 
of our computational technique we also consider the particular case 
of a Coulomb potential. Our result confirms the exp —C(logx)2 long- 
distance behavior first obtained by Schulz. 
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Models of one-dimensional (Id) interacting electrons are relevant as 
testing grounds for new ideas and can help us to understand systems in 
higher dimensions. In particular they are useful to describe the behavior of 
strongly anisotropic physical systems in condensed matter, such as organic 
conductors [p]] , charge transfer salts [@ and quantum wires || . Probably the 
two most widely studied Id systems are the Hubbard model j|] and the so- 
called "g-ology" model ||. They are known to display the Luttinger liquid 
H behavior characterized by spin-charge separation and by non-universal 
(interaction dependent) power-law correlation functions. In these models 
one usually considers short-range electron-electron interactions. This picture 
works well for conductors in which the screening between adjacent chains 
reduces the range of interactions within one chain |7|]. On the other hand, 
as the dimensionality of the system decreases charge screening effects are 
expected to become less important and the long-range interaction between 
electrons seems to play a central role in determining the properties of the 
system. This assertion seems to be confirmed by experiments in GaAs quan- 
tum wires and quasi-ld conductors ||. From a theoretical point of view 
the effects of long-range interactions has been also recently discussed in con- 
nection to a variety of problems such as the Fermi-edge singularity ||, the 
insulator-metal transition |I| and the role of the lattice through umklapp 



scattering and size dependent effects [O 



The effect of Coulomb forces on the single-particle Green function and 
on charge- density correlations in Id systems have been previously investi- 
gated in a pioneering work by Schulz W2\ using the conventional bosonization 



method [T3|]. The purpose of the present paper is to compute the electronic 
Green function in the presence of long-range interactions using an alternative 
path-integral bosonization technique previously developed in the context of 
quantum field theories We consider a non-local version of the Thirring 
model [15] described by the following (Euclidean) action 



S 



J d 2 x$0V - J d 2 xd 2 y [v m (x,y)Jo(x)Jo(y) + V {1) (x,y)J 1 (x)J 1 (y)] , (1) 



where x = (x, t x ), and Jo = ^70^ = P and J\ = \I>7i\I/ = j are the charge 
density and current density operators respectively. Let us also mention that, 
for simplicity, we shall take vp = 1 and pp = throughout this article. The 
functions V^(x,y) are forward-scattering potentials. Setting V( ) = V(i) = 



1 



Y^ 2 (x — y) one gets the covariant and local version of the Thirring model 
usually studied in the context of (1+1) QFT's. On the other hand, the choice 

V( )(x,y) = V(|x-y|) 5(r x -r y ), 

V { i){x,y) = 0, (2) 

yields the simplest version of the Tomonaga-Luttinger (TL) model with an 
instantaneous distance dependent potential and no current-current fluctu- 
ations. The main purpose of the present work is to study this case for a 
long-range electron-electron potential, different from the Coulomb interac- 
tion, of the form 



nix-y|) 



e 2 d~ 
|x-y 



1-e 



6(t x - T, 



(3) 



where d is a length scale and e is a small real number (0 < e < 1). This in- 
teraction has been previously analyzed by using renormalization group tech- 
niques in order to explore the possibility of having Luttinger liquid behavior 
in more than one dimension [[16] and in connection to the vacuum structure 



of the non local Thirring model [17 . 



Before considering the "super long range" potential described above, we 
shall sketch our functional approach to bosonization in a general case. As 
shown in ref. [JTSj the partition function of the model defined by eq. ([I]) 
can be solved, in the path-integral framework, by means of a chiral change 
of fermionic variables. This procedure allows to obtain a bosonized effective 
action in terms of scalar fields that are naturally identified with the collective 
modes (plasmons) of the system. Since this method has been described 
many times in the literature, here we shall skip the details (the interested 
reader is refer to [14}|, [15 and references therein). Let us only say that 
using a Hubbard- Stratonovich-like identity, which amounts to introducing an 
auxiliary vector field A^, the partition function of the general model given 
by ([I]) can be written in terms of a functional fermionic determinant as 

Z = JVA^ det{0 + gA) e~ s[A \ (4) 
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with 



S[A] = \J d 2 xd 2 yV { -^x,y)A fM (x)A fM (y), (5) 
where V^ 1 is such that 

f<PzVrf(z,x)V {lt) (y,z) = P(x-y). (6) 

Decomposing A^{x) in longitudinal and transverse pieces 

A^x) = e^d^x) + <9 M r?(x), (7) 

where and rj are boson fields (to be associated to the normal modes of the 
system) and applying, as anticipated, functional bosonization techniques to 
express the fermionic determinant in terms of <fi and i], one finally obtains 



Z = N 



J V(f)Vr] e~ Sb °% (8) 



where N is a normalization constant and Sb os , already written in momentum 
space, reads 

s bos = 7^-2 J d2 P [kp)k-p)Mp) + v(p)v(-p)B(p) + 4>(p)v(-p)C( P ) 



with 



and 



(9) 

A(P) = 1 -P 2 + \[V^{p)p\ + V a) \p)pH (10) 

m = l[v {0 ;(p)pi + v w l (p)pi], (ii) 

C(p) = [V (0) 1 (p)-V w 1 (p)}poPi, (12) 



A = C 2 (p) - AA(p)B(p). (13) 



Let us now focus our attention on the one-particle fermionic propagator 
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<.(*)*(,) >=( Gi » >v) G <f>). (14) 

The above depicted bosonization scheme can be easily employed to evalu- 
ate this expression. First of all it is straightforward to verify that the non-zero 
components of the Green function factorize as 



Gr, l (x, y) = G { r\(x, y)B RyL (x, y), (15) 



where G RL {z) = 2n (T^+z 2 ) ( Tz describe free right and left propagation, 

and 

B rA^v) = - r-4— . (16) 

J D^Drj e~ b bos 

In this equation one has defined 

Sr,l(x, y) = --^ J d 2 p [±4>(p) + ii)(p)] (e-*-* - e -*v). (17) 

At this point we notice that the evaluation of B R ^{x,y) can be carried 
out by means of a convenient shift of the fields <fi{p) and fj(p). This standard 
procedure gives 

Br,l(x, y) = exp[I R , L (x, y)\, (18) 
where I Rj l is a functional of the potentials given by 

IrA x ^v) = — 2 / rf2 P sin2 [|p- (^-y)] 

x |V(Q)(P)V(1)(P)P 2 + [V(o)(p) ~ V(i)(p)](PD T»Pi) 2 
fp 2 [V(i)(p)p§ + V(0)(p)p 2 ]+P 4 

This result gives a very general expression for the 2-point electronic corre- 
lator as functional of density-density (V(o)(p)) and current- current (Vm(p)) 
interaction potentials. It can be used as an alternative route to analyze the 
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effect of long-range interactions on a Luttinger system. Of course, if one con- 

'*• 2 

siders the local and covariant case V( )(p) = V(i)[p) = \ one easily obtains 
the well-known Thirring behavior: 

B T R h [ Trin \x -y)<x\x-y \-WM 2 /HWM . ( 20 ) 



Let us now proceed with our main task and specialize formula (|19D for the 
TL model defined by equation (|j). Having V^(p) = and V(o)(p) — V{o)(Pi) 
greatly simplifies the integrand in (|l9l) which allows to perform the integral 
in p . Moreover, since we are concerned with the equal-time correlator, we 
take t x = T y . In this case there is no distinction between I R and II and one 
gets 



Ir,l(\x - yU = r y ) = J(|x - y|) = -2 jT ^ 



_2 v /^(p) + l 2_ 
x[l-cos(|x-y|.p)] (21) 



where we have used pi = p. It is convenient to check the validity of this 
equation by considering a non trivial potential such as the Coulomb interac- 
tion 



e 2 



n|x-y|)= , |9 = (22) 
x — y\ 2 + a 2 



with Fourier transform given by 

V (p) = 2e 2 K (rfp). (23) 

where d is a very small length scale. This computation is interesting for 
two reasons. Firstly, one should stress that equation (|T9|) has been checked 
only for local (short-range) potentials: the covariant Thirring model (see 
eq. (|20D ) and the short range TL model fTSf . On the other hand, the electronic 
correlator for a Coulomb potential has been computed by Schulz [O using 



standard (operational) bosonization [|T^] and thus our computation will give 



an independent confirmation by means of a different approach. To achieve 
this goal we consider a large distance approximation of (|19|) that consists in 
inserting |x — y| _1 in the lower limit of (|2l|) and disregarding the cosine term, 



which is valid for -? >> , , , we obtain: 
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Gr,l(* ~ y) = ±iC(A)sign(x - y) exp 



3/2' 



-a log 



(3d 



(24) 



where a and (3 are two constants that depend on e 2 and C(A) depends on 
an ultraviolet cutoff A. This is our first non trivial contribution. We have 
obtained, as anticipated, the behavior previously found by Schulz fll2"| . In 
passing we note that a similar behavior has been previously observed by 
considering a Coulomb interaction between the core-hole and the conduction 



electrons 118 



Now we undertake the computation of G\r,z,(x— y) for the so called "super 
long-range" potential given by (|3|), whose Fourier transform is 



V(p) = 2e 2 cos -6 r(e)(d|p|) 



7T 



(25) 



Inserting this expression in fl2~l~|) one gets 



no 



|q| 



|q|" e + 2 



- 2 



2V|qh e + l 



[1 - cos(^|q| 



(26) 



where it has been natural to define the new variable 
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-e 2 cos 



7T 



7T 



r(e) 



i 1 ^ x-y 



(27) 



Note that we have also introduced a small quantity A to be set equal to zero 
at the end of our computation. Although I(£) is infrared convergent, it is 
convenient to split it into two pieces: /(£) = — = 0) with 



1 r 
^ = 2/ 



d|q| 



|q[- e + 2 
q|- e + 1 



cos(£|q|). 



(28) 



The first ^-independent piece can be very easily computed yielding 

I^Z = 0) = - V e/2 - log A - -(log 4 - 1) + 0(A e/2 ). 
e e 



(29) 
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The evaluation of ii(£) is more involved, as expected. It can be achieved 
by splitting the integration region into two sub-intervals (A, 1) and (l,oo), 
performing then series expansions of the corresponding square roots factors 
and finally integrating term by term. Putting all this together and taking 
the limit A — > one finally obtains 



+ ||e| e/2 r(-f ) cos(fe) + log |e| + (C + \ - \ log 4) (30) 



where §(a,b,z) is Kummer's confluent hypergeometric function, T(a,x) is 
the incomplete Gamma function, and Ci(z) is the cosine integral function. 
We have also defined the coefficients \i n — (n + |)e, v n — (n + l)e, and 

a n = ^ r (i)r(n+i) " F° rmu l a P0| ) is our main result. Indeed, by recalling that 

^(0 = ^(0 expJ(0, (31) 

one sees that fl30D gives a complicated but analytical and exact result for 
the fermionic propagator of the TL model in the presence of the long-range 
interaction given by equation (]3|). 

Of course, it is now interesting to study the long-distance behavior of our 
solution (|£| ^> 1). A careful evaluation of the dominant contributions to 
eq.(|30) for this case yields 



~r(-|) C0S (^)|^/ 2 



(32) 



Note that, exactly as it happens with Schulz's solution, the above function 
decays faster that any power-law. Moreover, the larger the range of the 
potential (larger e) the faster is the Green function decay. Let us also mention 
that the definition of the long-distance regime is different in both cases (see 



eq.(|27D), since the quantity e determines the range of the potential. This 
does not prevent us from comparing both long-distance decays for definite 
values of e. A simple numerical comparison allows to verify that the results 
are increasingly similar when distances increase. 
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It is also interesting to notice that the long-distance regime can be ob- 
tained following a much simpler route, analogous to the one employed to 
get this regime for the Coulomb case (see the paragraph preceding eq.(^|)). 
However, we must mention that the condition that allows to drop the os- 
cillating integrals leads, in the present case, to a result that is only valid 
for small e. Under this condition, replacing (|3|) in (|2lD and disregarding the 
cosine term gives 



Gr,l(0 k ±* sign(£)exp 



(33) 



which coincides with (|32|) for t <C 1. We want to stress here that this 
agreement is an important consistency check of our exact result since it has 



been derived in a very straightforward way, without using equation (31 ). 



In summary, we have shown how to compute, through path-integral meth- 
ods, the equal-time fermionic one-particle Green function in a simple version 
of the TL model. In particular we have obtained a rather involved but exact 
analytical expression for this propagator in the case of a " super long-range" 
potential (see (|3|)). This result is indeed of academic interest, since most of 
the explicit results found in the literature actually correspond to long dis- 
tance regimes (see for instance |12|). In passing, and as a consistency check 
of our computation, we have given a path-integral derivation of the Green 
function corresponding to the Coulomb well-known result previously 

found by Schulz. 
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